Search for Fermion Actions on Hyperdiamond Lattices 
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Fermions moving in a two-dimensional honeycomb lattice (graphene) have, at low energies, chiral 
symmetry. Generalizing this construction to four dimensions potentially provides fermions with 
chiral symmetry and only the minimal fermion doubling demanded by the Nielsen-Ninomiya no-go 
theorem. The practical usefulness of such fermions hinges on whether the action has a necessary 
set of discrete symmetries of the lattice. If this is the case, one avoids the generation of dimension 
three and four operators which require fine tuning. We construct hyperdiamond lattice actions with 
enough symmetries to exclude fine tuning; however, they produce multiple doublings. The limit 
where the actions exhibit minimal doubling does not possess the requisite symmetry. 

PACS numbers: 12.38.Gc 



I. INTRODUCTION 

When naive, chirally symmetric discretizations of spin 
one-half fermions are used, the famiUar phenomenon of 
doubhng arises: even though only one Dirac fermion is 
attached to each lattice node, the continuum hmit de- 
scribes several fermions. This is a problem for most ap- 
phcations, hke QCD, where only a small number of light 
quarks exist. For the naive fermion action in four space- 
time dimensions, doubling leads to 2 = 16 fermions. 
More complicated discretizations can reduce the number 
of doublings but not eliminate them completely, at least 
not while preserving exact chiral invariance. This is es- 
sentially the content of the celebrated Nielsen-Ninomiya 
"no-go" theorem [l|. 

In the past, Karsten [3] and Wilczek Q] proposed 
fermion actions with good chiral symmetry and the min- 
imal fermion doubling allowed by the no-go theorem. 
These actions were not successful in practice because they 
fail to have the correct continuum limit. The cuplrit of 
the problem is the lack of enough symmetry to forbid 
the existence of dimension three and four operators in 
the low-energy effective theory describing the lattice ar- 
tifacts. Thus, these actions require additional terms with 
coefficients fine tuned in order for the correct continuum 
limit to arise as the lattice spacing goes to zero. Usually 
fine tunings of this kind are very difficult to implement. 

Recently, Creutz 4] and, shortly after, Borigi 5] sug- 
gested a two-parameter class of actions with minimal 
doubling that generalized graphene to four space-time 
dimensions. Their action can be interpreted as being 
defined on a four-dimensional generalization of the two- 
dimensional honeycomb lattice (which is sometimes re- 
ferred to as a hyperdiamond lattice). It was subsequently 
pointed out [^] that this action lacks enough symmetry 
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to prevent dimension three and four operators from be- 
ing generated. The hyperdiamond lattice, however, has 
a symmetry under the permutation of any two out of five 
distinct axes. This 5*5 symmetry (or a large enough sub- 
group, such as the group of cyclic permutations Z5), if 
possessed also by the action, would eliminate all dimen- 
sion three and four operators that cause the fine tuning 
problem. A natural question then is whether the Z5 sym- 
metry can be recovered for special values of the action 
parameters, as suggested in Q, or whether, failing that, 
a modified Z5 symmetric action can be engineered. In 
order to clarify this issue, we construct a Z5 symmetric 
action on a hyperdiamond lattice. This action exhibits 
multiple doublings, and attempts to eliminate these dou- 
blings lead to other undesirable features. 



II. HYPERDIAMOND LATTICE ACTION 

A. Hyperdiamond lattice 

We start the construction of the four-dimensional hy- 
perdiamond lattice by defining a set of five distinct vec- 
tors {e"} satisfying 



e^-e'^ = 



1, for a — (3 
cos 9, for a ^ (3 



(1) 



Straightforward linear algebra shows that the mutual an- 
gle between these vectors is given by cos 9 — —1/4. Con- 
sequently we have the relation J2a=i ®" = 0. We can 
take the vectors e" to be 

el = -( V5, V5, VS, 1), 



= i( V5,-%/5,-V5, 1), 

= i(-V5,-V5, Vb, 1), 

= i(-V5, \/5,-V5, 1), 

= -( 0, 0, 0, 1). 
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The first four elements of this set form a basis whose dual 
basis {e^} is given, up to normalization, by 
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-e^ 


62 = 


9 

= e 


-e^ 


63 = 


= e^- 


-e^ 


64 = 
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(3) 

The hyperdianiond lattice is formed by two sublattices. 
One (the "L-nodes" ) are generated by integer combina- 
tions of e^, that is, the points x^ = X)u=i ^'^^M' with 
integer x'^ . The second sublattice (the "R-nodes" ) is ob- 
tained from the first by translating it by e^, so that it is 
formed by the points x^ — J2u=i ^^^t^ + 6^. Notice that 
the four coordinates x'^ do not specify a node uniquely: 
we need to specify, in addition, whether the point is an 
L-node or an R-node. 

The five nearest neighbors of an L-node with coordi- 
nates x^ arc the R-nodes with coordinates 

I it/ I -L-ii ^ Jj . lij J • {•Jj ^ Jj I _L , it/ , Jj j . 

{x^,x^,x^ + l,x'^), {x^,x^,x^,x'^ + 1), and 
{x\x\x\x^), 

as is easily verified. Similarly, the nearest neighbors of an 
R-node with coordinates [x^ , x"^ , x^ , x'^) are the L-nodes 
with coordinates 
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The separation between one node and its nearest neigh- 
bors are the five vectors e". The symmetric distribution 
of these vectors means that the nearest neighbors are also 
distributed symmetrically around each node. In fact, not 
only the nearest neighbors but the whole hyperdiamond 
lattice is symmetrical under any permutation of the five 
e" axes. Hence the hyperdiamond lattice has an 6*5 sym- 
metry. These permutations are of crucial importance, 
and so we examine them carefully. The permutations 
involving only the first four axes are described in the co- 
ordinate system used here by the same permutation of 
the x^ coordinates. For instance, exchange of the e^ and 
e^ axes is the mapping 



{x\x\x\x^ 



] ' ' ' X ~~ let/ ,it/ • tLi ,tt/ ] ^ 



(4) 



for both L-nodes and R-nodes. An exchange involving 
the fifth axis is slightly more complex. The e^ and e^ 
exchange, for instance, takes the L-node 

{x^,x^,x^,x'^) — > x^ {-[x^ +x^ + x^ +x'^),x'^ ,x^ ,x'^), 

(5) 



but takes the R-node 

{x^,x^,x^,x'^) — > x^ {-{x^+x^+x^+x'^) + 



l,x\x\x^) 



(6) 



B. Hyperdiamond action 

We can build a free fermion action on the hyperdia- 
mond lattice by attaching a two-component left-handed 
spinor 4> and right-handed spinor (f) to each L-node, and a 
right-handed spinor x and left-handed spinor x to every 
R-nodc. The action 



^ = E 



M=l 



.(^■e'^Xx 



Xx+^i CT • e^ 0^) 



i^x^-^^Xx-Xx^-^^' 



(7) 



with CT — (a, — i), a — ((?,«), describes fermions hopping 
to nearest-neighbor sites with equal probability in all five 
directions. Since (px and x^ live on L- and R-nodes, re- 
spectively, the last two terms describe the hopping along 
the e^ direction even though the coordinate x is the same 
for both fields.^ 

The action in Eq. ^ is not invariant under the full S'5 
symmetry group of the lattice. However, it is invariant 
under the subgroup of transformations which can be built 
by an even number of elementary permutations of two 
axes, known as the alternating group, A^. Under an A5 
transformation, the fields transform as 
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where the permutation operators are 
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and the coordinates transform as in Eqs.fll])-®. For ex- 
ample, under a transformation involving a permutation 
of the axes e-'^ and e^, as well as a permutation of e-'^ and 
e'^, the field 6,- transforms as 



P12P 



12 J^^13 'Px3,Xi,X2,X4 



(10) 

The transformation of the action in Eq. ([7]) under per- 
mutations of the five axes follows from 



7^_„/3 



Pap (J ■ e" Pap = cr • e'^ 

PaP (J ■ e" PaP = CT 

Pap a ■ e^ PaP = (J ■e' , a,ii 4^1 

Pap o ■ e'' Pap = a ■e' , a,ii ^-^ 



(11) 



^ In momentum space, this aetion reads 






(^sin(p^)e'" -7- (^cos(p^)e'' + e^ j -775 
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with 



'/'P=(/). V'p=(0p,Xpj. and 7m = U 



cr. 



which, in turn. foUows from 



Dirac spinor ■0: 



2a-bb-a 



{a-ab-a+b-aa-a)b-a = a-a+b-aa-ab-a, (12) g^^, — _ V^ 



valid for any unit vectors a and b. Being careful to re- 
member that L-nodes and R-nodes transform differently, 
and that 0, (j) sit on L-nodes, while x, x sit on R-nodes, 
one can verify that the action, Eq. ([7]), is invariant under 
the subgroup of even permutations of S5 , the alternating 
group, ^5. 

There are only two chirally symmetric relevant op- 
erators in the continuum that are ^5 symmetric: 
*Z]aV'(a;)e" • 7'0(a;), and J^a'^i^)^"' ' 775tp{x). Both 
of these operators, however, identically vanish because 
the five basis vectors sum to zero. In fact, we need not 
require invariance under the full A^ symmetry to ex- 
clude such operators, the cyclic permutation subgroup 
Z5 C A5 is the minimal subgroup of S5 that excludes 
relevant operators. While the action exhibits enough 
symmetry for a good continuum limit, it does not, how- 
ever, have minimal doubling. In fact, it has, in addi- 
tion to the pole at p^ = 0, where the action reduces 
to the Dirac form, several others poles at finite p^, e.g. 
Pi — — p2 — —P3 — Pi = cos~"'^(— 2/3). Let us now 
consider actions with minimal doubling and argue that, 
unfortunately, they seem to lack the crucial Z5 symmetry. 



III. COMPARISON WITH OTHER ACTIONS 

Creutz m, and later Boriqi [^, suggested minimally 
doubled fermion actions based on similar considerations 
to the ones discussed above. In fact, the Borigi-Creutz 
action can be written in the form 



'sc 



^E 



^(V;._^e^.rv,-^,+^e^-rt^,) 



-H^.e^.rVx-V'xeS.rtV', 






+ X,E-e5 0,~(/.,S-e^X. 



(15) 



, (13) 



with S = (tr, —1) and E = (tr, 1). 

Despite the similarity between S in Eq. ([7]) and Sbc in 
Eq. P^ . there are important differences. First of all, for 
generic B and C, the vectors e" are not symmetrically 
arranged and thus break the Z5 symmetry of the lattice. 
Only the case B — l/\/5 and C — 1 satisfies Eq. ([T]), up 
to the overall normalization. A further difference is that 
there is an extra factor of i in the fourth components of 
S, E compared to cr, ct. The reason both actions are able 
to reproduce the Dirac propagator despite this additional 
i factor is that while S generates a pole at p^ — 0, the 
poles of Sbc are at finite p^, where it duly reduces to 
the Dirac action. The different pole positions also ex- 
plain why the position of S, E appear shifted from the 
positions of ct, ct, even though both actions reduce to the 
Dirac action around their poles. The differences in spin 
structure do not spoil the spinor part of the 6*5 transfor- 
mations provided that B = l/-\/5j and C = 1. Finally, 
Sbc contains twice as many terms as S: the additional 
terms comprise the last two lines of Eq. (IT5|) . These extra 
terms can be visualized as non-nearest neighbor interac- 
tions in the hyperdiamond lattice. To see this, take the 
case where x — (0, 0, 0, 0). As shown in Fig. [U the term 
Xx-i '^' ^^ 4>x, for instance, describes a hopping from the 
point X = to — ei + e^. Most importantly, the addi- 
tional terms in the last two lines of Eq. P^ break the Z5 
symmetry under coordinate transformations. Without 
the additional symmetry, relevant and marginal opera- 
tors will be generated by gauge interactions Q. 



IV. MODIFIED ACTIONS 



where F^ = (7,^74) and the vectors e" are defined in 
terms of two parameters B and C as 
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In order to compare this action to Eq. ([7]) , let us write it 
in terms of the left- and right-handed components of the 



Guided by our construction of an A^ symmetric action 
on the hyperdiamond lattice, we can try to answer the 
question: is there a way of combining the good features 
of the Borigi-Creutz action (minimal doubling and chiral 
symmetry) while avoiding its pitfall (lack of Z5 symme- 
try)? We have been unable to find such an action which 
answers this question. For instance, one possibility is to 
demand invariance under Z5 and hence drop the terms 
in the last two lines of Eq. P^ . This corresponds to in- 
cluding chiral projection operators in Eq. p3p . We must 
further choose B = l/^/b and C = 1. The resulting ac- 
tion is Ac, symmetric with a pole at p^ = 0, and several 
other values of p^. The pole at p^ = 0, due to the mis- 




FIG. 1: The two-dimensional analogue of the hyperdiamond 
lattice. The long (blue) lines denote the dual basis, {ei,e2}, 
as defined in Eq. ((3)|. In both figures, the (black) circle sits on 
an L-node, and denots the field <^o,o. The (red) squares denote 
various x fields sitting on R-nodes. Top figure: the three 
nearest neighbor fields, Xo O' Xi o ^'^'^ Xo 1 1 connected by thick 
(black) links. Lower figure: the three (non)-nearest neighbor 
fields, i^Q g, x_i ^iid Xo -II connected by thick (black) links, 
present in the Sbc action. 



the Dirac equation in the naive continuum hmit. Fermion 
poles of this form, sometimes called mutilated fermions, 
also appear in other attempts at formulating fermions on 
non-hypercubic lattices [TI, '8| . 



V. CONCLUSION 

By analogy with fermions on a honeycomb lattice, we 
construct a hyperdiamond fermion lattice action. This 
action possesses a large subgroup (^5) of the 6*5 sym- 
metry of the hyperdiamond lattice, and consequently 
enough symmetry to avoid fine-tunings in taking the con- 
tinuum limit. The hyperdiamond action has, however, 
more than the minimal amount of fermion doubling. 

We investigate the lattice actions proposed by Creutz 
and Borigi, and find that, while similar to the hyperdia- 
mond construction, these actions do not have at least a 
Z5 symmetry. Hence the continuum limit of these theo- 
ries will require fine tuning. Modification of these actions 
to enforce Z5 symmetry from the outset leads to unde- 
sirable effects such as additional poles, and a Lorentz 
non-symmetric continuum limit. 

The goal of our search is to produce a minimally dou- 
bled action with a good continuum limit. We find, how- 
ever, an intricate balance needed to break hypercubic, 
parity, and time-reversal symmetry in order to obtain 
minimal doubling, while, at the same time, preserve or 
invent additional symmetries necessary to avoid fine tun- 
ings. The requirement of Z5 symmetry on a hyperdia- 
mond appears incommensurate with minimal doubling. 
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placed i factors and S matrices relative to Eq. ([7]), has 
the wrong form. The action around p^, ~ behaves like 

■0^,(^7 • fc -I- 7574fc4)i/'/c, for suitably defined physical par- 
ticle momenta k^^ = kp,{pu)- This form fails to reproduce 
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